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Overview 



• Motivation and applications 

• Multiscale materials modeling 

• Atomic structure 

• Interatomic potentials 

• Simulations of lattice thermal conductivity 
for ZrB 2 and HfB 2 

• Comparison to experiments 
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UHTC for Sharp Leading Edges 



Sharp leading edge for hypersonic aircraft 

- Enhances vehicle performance 

- Improves safety 

Higher temperature requirements 

- Shuttle RCC leading edge: T~1650C 

- Sharp leading edged vehicles: T>2000C 
UHTC advantages for sharp leading edges 

-Good mechanical properties 

- Oxidation resistance 

- High themngl conductivity 

- Effective thermal radiation 

- Thermal shock resistance 



Leading edges of hypersonic vehicle 



SEM image of UHTC microstructure 


Gasch, Johnson and Marschall, J. Am. Ceram. Soc. 91, 1423, (2008) 
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Multiscale Modeling of UHTC 



• Framework integrates three methods 

• Multiscaie framework for ZrB 2 and HfB 2 : 

- Ab initio - fundamental chemistry, 
electronic structure impact on basic 
material properties 

- Atomistic - thermal/mechanical 
properties, adhesion and thermal 
resistance of grain boundaries, fracture 

- Continuujv- macro properties, 
thermal/mechanical analysis of 
microstructure 

• This talk focuses on atomistic methods 

- Development of interatomic potentials 

- Lattice thermal conductivity simulations 

- Other topics presented elsewhere 



JL, Daw, Squire and Bauschlicher, (2012), submitted 
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NASA 




Atomic Structure: ZrB 2 and HfB 2 









Alternating layers of 
Zr/Hf (red) and Boron (gray) 


r*Y#Y*Y«r*Y* 

r^T^rir^T^ 


Graphitic Boron layers 
with Zr/Hf over each ring 
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Fundamental Properties: ZrB 2 & HfB 2 


Electronic Spectra 



E-Ef(eV) 

Electronic properties 
essentially identical 


Vibrational Spectra 



Frequency (cm-1) 


Acoustic modes carry heat . 
Optical modes are resistive . 


JL, Bauschlicher and Daw, J. Am. Ceram. Soc., (2011) 
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Tersoff Bond Order Potential 


• Two body terms (A,X, B, jn) energy 

/•- Zl// 

i*j 

f^(d) = A ij e^p(-A iJ d) 
fA ] (d) = -B ij exp(-v ij d ) 


• Bond order (P, ^ 3 , n, m) 

fy=c :i+A"'Q)~ /2ni 

iij = Z fc Jl ( r ik ^YijkSi Wqk ) eX P [^ 3 , (djj ~ d ik )”' ] 

k*i,j 

• Angular function (c, d, h) 

g t (0)= 1 + cf / df - cf /[df + (h t - COS 0 2 ] 



Interatomic Energy 


Daw, JL and Bauschlicher, Comp. Mat. Sci., (2011) 
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NASA 
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First Step: Zr Potential 


• Zr potential exists 

• Developed new Zr potential 

• Fit to ab initio database of crystal structures 


Property(units) 

Target 

New 

WM2 

a 0 (FCC) (A) 

4.530 

4.510 

4.532 

E 0 (FCC) (eV) 

-6.160 

-6.159 

-6.127 

B(FCC) (eV/A 3 ) 

0.578 

0.5899 

0.6011 

B’(FCC)(eV/A 4 ) 

-0.8160 

-1.635 

-1.948 

C 11 (FCC)(eV/A 3 ) 

0.7740 

0.6885 

0.7404 

C 12 (FCC)(eV/A 3 ) 

0.4810 

0.5405 

0.5314 

C 44 (FCC)(eV/A 3 ) 

0.3560 

0.5307 

1.395 

E vac (FCC)(eV) 

2.500 

6.072 

8.338 

a 0 (HCP) (A) 

3.230 

3.159 

3.231 

E 0 (HCP) (eV) 

-6.180 

-6.242 

-5.826 

E 0 (BCC) (eV) 

-6.050 

-6.159 

-5.960 



FCC 


HCP 


BCC 


Williame and Massobrio, PRB 43 (1991), 11653 
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Second Step: Boron Potential 


• No pubished Boron potentials 
•Boron is electron " deficient " 

• Boron may be "frustrated" 

• Fit to simple structures 


Structure 

Property 

Target 

Fit 

Hex sheet 

a o 

2.91 

2.89 


E 0 

-5.15 

-5.08 


E 0 ” 

11.35 

7.98 

Tri sheet 

a o 

1.70 

1.81 


E 0 

-5.71 

-5.75 


F ” 
c 0 

21.73 

27.06 

SC 

a o 

1.88 

1.84 


E 0 

-5.33 

-5.21 


E 0 ” 

24.50 

24.51 

FCC 

a o 

2.86 

2.84 


E 0 

-5.07 

-5.22 


E 0 ” 

21.85 

12.28 



P-Boron (N=1 05) a-Boron (N=12) 



Crumpled Sheet 


Ogitsu et al, JACS 131 (2009) 1903 
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Third Step: ZrB 2 


• Zr-Zr parameters fixed 

• B-B parameters fixed 
•Zr-B fit to small database 

• Pot A = "new Zr" + B 

• Pot B = WM2 + B 

• Will Boron planes stay flat? 


Fitting Results 


Property 

Target 

Pot A 

Pot B 

a o(A) 

3.170 

3.143 

3.140 

C 0 (A) 

3.550 

3.547 

3.547 

E 0 (eV) 

-21.70 

-21.29 

-21.55 


Potentials 
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Stable, multilayered system 
with flat, hexagonal Boron sheets! 


O Boron O Zirconium 





Lattice Thermal Conductivity 




Green-Kubo thermal conductivity tensor 



Heat current J(x ( -,V/1 energy e lf stress-tensor S t 


7 = 1 
V 

7 = 1 

V 


Z e < v < -Zv, 


Z e ^ + ^Z(4-( V . +V ;))- 


X s 


JL, Daw and Bauschlicher, J. App. Phys, (2011) 
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Heat Current Correlation Function 
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• Monoatomic systems (e.g. Si) have monoatomic decay 

• ZrB 2 has longer period than HfB 2 at T=300K 

• ZrB 2 at T=1000K has longer period than T=300K 
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Power (ar bitary units) 


Correlation Function Power Spectra 



• Correlations oscillates with metal-B optical modes 

• C xx and C yy oscillate with in-plane mode frequency 

• C zz oscillates with out-of-plane mode frequency 


K(W/(m.K)) 



Lattice Thermal Conductivity: ZrB 2 



8 independent, 10 ns simulations, T=300K 
8x8x16 unit cell, 12,255 atoms 
k xx =60 W/(m.K), k zz =40 W/(m.K) 
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K(W/(m.K)) 


Lattice Thermal Conductivity: HfB 2 




• 8 independent, 10 ns simulations 

• 8x8x16 unit cell (12 atoms)= 12,255 atoms 

• k x =76 W/(m.K), k z =65 W/(m.K) 
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k (W/(m.K)) 



Thermal Conductivity vs Temperature 



•8 independent, 10 ns simulations for each point 

• Data fit to 1/T curves 
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Experimental Data Comparison 



• Polycrystalline ZrB 2 

• K e =33 W/mK, K, at = 22 W/mK 

• K lat ~ 0-3K tot 

• Single crystal ZrB 2 

• k xx = 140 W/mK, k zz = 100 W/mK 

• 1 sample, 1 measurement 

• defects uncharacterized 

• k xx = 45 W/mK , k zz = 30 W/mK 

• Data needed for ZrB 2 and HfB 2 

• Simulation data reasonable at 300K but 

too low for higher T 


Polycrystalline ZrB 
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Zimmermann, Hilmas, Fahrenholtz Dinwiddie, Porter, Wang , J. Am. Ceram. Soc., (2008) 
Kinoshita, Otani, Kamiyama, Amano, Akasaki, Suda, Matsunami, Japan. J. App. Phys., (2001) 
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Conclusions 


• Atomistic simulations for ZrB 2 and HfB 2 : 

• Developed first interatomic potentials for UTHC 

• Lattice thermal conductivity using Green-Kubo formalism 

• Heat current correlation function oscillations 

• Thermal conductivity versus temperature 

• Reasonable agreement with experiment 

• Modeling unanswered questions: 

• Interatomic potential fidelity 

• Lattice TC without potentials ( ab initio , Boltzmann,. . .) 

• Conducting versus resistive vibrational modes 

• Isotope and defect effects 

• Interface thermal resistance: grain boundaries * 

• Experimental unanswered questions: 

• Single crystal characterization and thermal conductivity 

• Electronic versus lattice thermal conductivity 


* JL, Daw, Squire and Bauschlicher, (2012), submitted 



Extra Slides 
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ZrB 2 Potential Curves 


Test Results 


Properties 

Ab Initio 

Pot A 

PotB 

Cn 

556 

365 

422 

C i2 

57 

156 

156 

^13 

113 

173 

171 

CO 

CO 

o 

419 

307 

320 

C 44 

234 

106 

119 

B 

233 

227 

240 

G 

226 

98 

118 

A(=C 33 /Cii) 

0.75 

0.84 

0.76 


Properties not included in fit 
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Lattice Constant (A) 


VASP = ab initio code 





Physics of Lattice Thermal Conductivity 


•K=pc v/ mfp 

• scattering restricts / mfp 

• Region I: k ~ T 3 

• dilute phonons 

• boundary scattering 

• quantum statistics 
•Region II: K niax 

• Region III: k ~ 1/T 

• high phonon density 

• phonon, pt. defect scattering 

• Region IV: K min , / mfp = “a” 



Clark, Surf, and Coat. Tech., (2003) 
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Summary 


• No atomistic simulations for ZrB 2 due to lack interatomic potentials 

• Potentials are prerequisite for atomistic simulations of mechanical and 
thermal properties 

• We developed such potentials for ZrB 2 

• ZrB 2 potentials give stable structures with flat, hexagonal B planes 

• We performed the first atomistic simulations for these materials 

• Lattice thermal conductivity was evaluated for single crystals 

• Reasonable agreement with experiments 

• Future/current work: 

• Grain boundaries: energetics and thermal interface resistance 

• Integration into multiscale framework 

• Potentials and applications for H/and HfB 2 


JL, Daw, Squire and Bauschlicher, (2012), submitted 
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